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ABSTRACT 

Rossi X-ray Timing Explorer (RXTE) All-Sky Monitor light curves of Circinus X-l 
show that intensity dips below the ~1 Crab baseline occur near phase zero of the 16.55-d 
cycle of the source. RXTE Proportional Counter Array observations carried out between 
1996 September 20-22 provided 60% observing efficiency for 48 hours around phase zero. 
These observations showed significant dipping activity during much of those two days. 
The dramatic spectral evolution associated with the dips is well-fit by variable and at 
times heavy absorption (Nh > 10 24 cm -2 ) of a bright component, plus an underlying 
faint component which is not attenuated by the variable column and whose flux is about 
10% of that of the unabsorbed bright component. A prominent iron emission line at 
6.4-6.6 keV is evident during the dips. The absolute line flux outside the dips is similar 
to that during the dips, with equivalent width increasing from as low as 44 eV outside 
dips to more than 400 eV during dips, indicating that the line is associated with the faint 
component. These results are consistent with a model in which the bright component 
is radiation received directly from a compact source while the faint component may be 
attributed to scattered radiation. Our results are also generally consistent with those of 
Brandt et al., who found that a partial-covering model could explain spectra obtained 
in ASCA observations of a low-to-high transition in Cir X-l. The relative brightness of 
the two components in our model requires a column density of ~2xl0 23 cm -2 if the faint 
component is due to Thomson scattering in material that mostly surrounds the source. 
We find that illumination of such a scattering cloud by the observed direct component 
would produce an iron Ka fluorescence flux that is in rough agreement with the flux 
of the observed emission line. We also conclude that if the scattering medium is not 
highly ionized, our line of sight to the compact source does not pass through such a 
scattering cloud. Finally, we discuss simple pictures of the absorbers responsible for the 
dips themselves. 
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1. Introduction 

Circinus X-l exhibits complex X-ra y variability 
during its 16.55-d intensity cycle (e.g 



Kaluzienski 



t al 1976| ; Dower, Bradt, k Mo rgan 1982; [Tcnnant 
1987|, 1988; [Stewart et al. 199l|, Shirey et al. 1996, 



1998, hereafter Papers I fc II; Shirey, Bradt, & Levine 

1999, hereafter Paper III ). Enhanced mass transfer 
near periastron of a highly eccentric binary orbit may 
be responsible for both flaring and dipping behavior 
flMurdin et al. 1980|; |Oostcrbroek et al. 1995|; |Brandt 
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i\. 1996 )r — Three type I X - ray burata accn from 



Cir X 1 during an EXOSAT obacrvation demonstrate 

that the system contains a neutron star with a weak 
magnetic field (Tennant, Fabian, & Shafer 1986). 

Rossi X-ray Timing Explorer (RXTE) All-Sky Mon- 
itor (ASM; 2-12 keV) light curves (see Papers I & II 
and [Shirey 1998 ) show that Cir X-l has maintained a 
bright baseline level of ~1 Crab (1060 [iJy at 5.2 keV) 
since early 1996. Dramatic flares generally begin dur- 
ing the day following phase zero, as given by the radio 
ephemeris (Stewart et al. 1991), of each 16.55-d cycle 
of the source. The flaring reaches as high as 3.7 Crab 
and typically lasts 2-5 days, but sometimes persists 
for most of the cycle. This behavior is similar to that 



observed in folded Ginga ASM data (Tsunemi et al 



1989) 



The RXTE ASM light curves of Cir X-l also show 
intensity points below the 1 Crab baseline near phase 
ze ro of about h alf of the cycles observed (see Fig. 3-1 
of Shirey 1998 ). Because ASM coverage is incomplete 
(about ten to twelve 90-s exposures per day), it is pos- 
sible that dips, lasting seconds to hours in duration, 
might actually occur in all cycles. During such dips, 
the X-ray intensity reaches as low as —80 mCrab. The 
dips most often occur during the half day before and 
after phase zero, but sometimes occur up to two days 
after phase zero. Occasionally, dips also occur at later 
phases in the ASM light curve, e.g., near day 4 of the 
cycle. 

X-ray dips are common in the light curves of X-ray 
binaries. In some cases, dips may be due to a true de- 
crease in the X-ray emission, e.g., due to an accretion 
disk instability in which the inner part of the disk 
is ejected (Greiner, Morgan, fc Remillard 1996; Bel- 



[oni et al. 1997 ). However, in most cases dips are at- 
tributed to absorption by matter passing through the 



line of sight to the X-ray emitting region (see White 
Nagase, & Parmar 1995 for a recent review). Deep, 
erratic dips indicate a compact emission region and 



structure in the absorber, probably near the outer 
edge of the accretion disk. Shallower, smooth dips 
occur in systems where the compact source is always 
obscured by the disk (due to a high inclination) so 
that only partial covering of an extended corona can 
occur. In some systems, X-ray eclipses by the com- 
panion star occur and are total or partial depending 
on whether the non-eclipse emission is dominated by a 
compact source or an extended region such a corona. 

An ASCA observation of Cir X-l at phase zero 
showed a low-to-high intensity transition (Brandt 



et al. 1996), apparently the egress from a dip simi- 
lar to those seen with the ASM. Brandt et al. used 
a partial-covering model to fit the spectra before and 
after the transition. A strong iron K edge in the low- 
state spectrum indicated that the obscuring matter 
had a very high column density, i.e., Njj > 10 24 cm~ 2 . 
Those authors suggested that scattering by electrons 
in circumstellar material might be the source of the 
faint (dip) spectrum. 

In Papers I, II, & III we presented RXTE results 
on the spectral and timing evolution of Cir X-l during 
quiescent and active phases, exclusive of the behavior 
during dips. Those results demonstrate that Cir X- 
1 at least sometimes exhibits Z-source behavior, in- 
cluding horizontal branch oscillations (HBOs) at un- 
usually low frequencies (1.3-35 Hz), normal branch 
oscillations (NBOs) at ~4 Hz, and only very low fre- 
quency noise (VLFN) on the flaring branch. In this 
fourth paper we focus on spectral evolution during 
dips observed with the RXTE Proportional Counter 
Array (PC A). 

2. Observations 

Observations of Cir X-l with the PCA were made 
during 1996 September 20-22, covering 48-hours around 
phase zero with —60% observing efficiency (see Fig. Q). 
These observations showed significant dipping for 
about half of the observed time and provided the op- 
portunity to study dip behavior in detail. 

The observations began about a half day before 
phase zero. The intensity of Cir X-l was very low (an 
extended dip) during most of the first half day of the 
observations, and then returned to the much higher 
"quiescent" level shortly after phase zero. Shorter in- 
termittent dips occurred during the —0.3 days follow- 
ing phase zero. The strongest dips all reached similar 
minimum levels of about 200-350 counts s _1 PCU -1 
(2-18 keV) after background (-10 c s" 1 PCU" 1 ) was 
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subtracted, where 1 Crab « 2600 c s" 1 PCU -1 . A 
second episode of significant extended dipping oc- 
curred from about 0.5 to 0.9 days after phase zero, 
followed by "active-state" behavior for the remain- 
ing half day of the observation. The active state is 
characterized by (1) a gradual increase of the aver- 
age intensity in the lowest-energy PCA band (which 
dominates the total PCA count rate) accompanied by 
a significant drop in the high-energy bands, and (2) 
individual flares on shorter timescales (less than a few 
hours) which show increased intensity in all PCA en- 
ergy bands. The former behavior results in gradual 
shifting of spectral tracks in color-color and hardness- 
intensity diagrams (CDs/HlDs) and is related to the 
16.55-d intensity cycle of Cir X-l. The latter behavior 
corresponds to motion along branches of the "Z" track 
in CDs and HIDs. For example, based on tracks in 
hardness-intensity diagrams and on timing properties, 
we find that the burst-like events near day 348.0 in 
Figure |l| correspond to motion up the flaring branch. 
The final two segments of Figure |l| (after day 348.1) 
show motion from the NB/FB apex up the normal 
branch to higher total intensity, and back down. The 
spectral and timing evolution during the active state 
of Cir X-l was discussed in detail in Papers II & III. 

3. Evolution of Hardness Ratios During Dips 

Three dips, each showing a complete transition 
from the baseline intensity level down to a deep min- 
imum, have been selected for further study and are 
labeled Dips 1-3 in Figure |l|. An expanded 2400-s 
segment of Figure [j] is shown for each of these dips 
in Figures ||, ^, and ||. The light curves in multi- 
ple PCA energy channels show that the transitions 
into and out of dips occur more rapidly at low energy 
that in higher energy bands; this is most pronounced 
for Dip 3 in Figure The hardness ratios thus ini- 
tially increase (harden) as the denominator falls more 
quickly than the numerator. Low-energy count rates 
are also the first to reach a minimum level well above 
background. This causes the hardness ratios to de- 
crease (soften) since the denominator approaches a 
constant while the numerator continues to decrease. 

The intensity and hardness ratios from Dip 3 have 
been used to produce the color-color and hardness- 
intensity diagrams (CD & HID) shown in Figure §. 
The spectral evolution during the dip produces tracks 
with dramatic bends in both the CD and the HID. A 
similar CD and HID (Fig. |^) has also been produced 



for the first three segments of Figure [l] (day 346.31- 
346.52), during which the intensity was mostly in an 
extended low/dip state but briefly reached higher in- 
tensity. 

In the HIDs of Figures || and ||, dips produce mo- 
tion to the left (toward lower intensity) and initially 
upward (harder), but turn downward (softer) as the 
count rate in the denominator approaches its bottom 
level. In the CDs where, by definition, two hardness 
ratios (hard and soft color) arc employed, motion is 
initially to the right and upward (harder in both col- 
ors but mostly in the lower channels) . When the low- 
est channel reaches the bottom level, the track turns 
to the left but continues upward. When intermediate 
channels stop decreasing and only the highest chan- 
nels still drop, the track finally turns downward. The 
intensity of the second and third lowest energy bands 
rapidly drop by a factor of two at time 1900 s in Fig- 
ure f|, just before settling to their lowest levels; this 
resulted in a gap in the the middle CD branch in Fig- 
ure |j| 

Evolution of the position in CDs and HIDs apart 
from dips is interpreted in Papers II & III as motion 
around a "Z" track as well as shifts of the entire "Z" 
pattern. The exact position of a dip track in the dia- 
grams depends on its starting point, i.e., the baseline 
intensity and spectrum. Although starting from dif- 
ferent positions, all dips below the 1-Crab baseline in 
PCA light curves of Cir X-l produce the same general 
shape in CDs and HIDs. However, the tracks may not 
be complete for weak dips. 

4. Evolution of Energy Spectra During Dips 
4.1. Spectral Model 

Cold (neutral) matter most effectively absorbs pho- 
tons at lower energies, so gradually increasing obscu- 
ration by cold matter produces dips that are more 
gradual at high energy than at low energy, as we ob- 
serve in Cir X-l. However the fact that strong dips in 
Cir X-l show a minimum intensity well above back- 
ground suggests the presence of a faint component 
unaffected by the dips. Thus, in modeling spectra 
of dips in Cir X-l, at least two components are nec- 
essary: (1) a bright component modified by variable 
and at times very heavy absorption and (2) a faint 
component seen through a constant and relatively low 
column density. From a study of the energy spectrum 
during non-dip observations (Paper III), wc found 
that the 2.5-25 keV spectrum of Cir X-l was generally 
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fit well with a model consisting of a multi-temperature 
"disk blackbody" plus a higher-temperature (~2 keV) 
blackbody. We find that a model based on these com- 
ponents also adequately fits the current spectra, so 
a disk blackbody plus blackbody is adopted for the 
bright component during dips. For simplicity, we as- 
sume that the unabsorbed spectral shape of the faint 
component is similar to that of the bright component, 
as might be the case if it is due to scattering, e.g., by 
a corona or stellar wind, of radiation coming directly 
from a compact central source back into the line of 
sight from other initial emission directions. Thus for 
the faint component we use the same model and pa- 
rameters as the bright component, but normalized by 
an energy-independent multiplicative factor / which 
is less than unity. We discuss the validity of this as- 
sumption below. 

The absorption column density for the faint com- 
ponent is assumed to remain constant throughout a 
dip, but the column density for the bright compo- 
nent is allowed to vary to produce the dips. The light 
curves during strong dips show a significant reduc- 
tion in intensity even at high energy (13-18 keV), 
indicating very high column densities (Nh > 10 24 
cm -2 ). Photoelectric absorption is the dominant pro- 
cess responsible for dips at low energy, but the pho- 
toelectric cross-section decreases with energy. Thus, 
above about 10 keV, Thomson scattering (with an 
approximately energy-independent cross section) be- 
comes dominant. Both photoelectric and Thomson- 
scattering factors are included in the attenuation cal- 
culation for the bright component. For simplicity, the 
effective hydrogen column density and the Thomson- 
scattering (electron) column density are taken to be 
equal; the actual value of the latter will be slightly 
higher (by ^20%) due to atoms (mainly helium) with 
more than one electron. The complete model used 
can then be expressed as: 



F = 



CX P (-cr p/l iVH (1) )cxp(-(TThiVH (1) ) 



M - 



exp 



(~a ph N H {2) )f} M, (1) 



where F is the observed flux, <j p h and gtk are the 
photoelectric and Thomson cross-sections, Nh and 
Nh^ 1 are the effective hydrogen column densities of 
the bright and faint components respectively, / is the 
ratio of the unabsorbed flux of the faint component 
to the unabsorbed flux of the bright component, and 
M is the disk blackbody plus blackbody model. For 
all model components we assume solar abundances 



( Anders fc Ebihara 1982j ) , and we used the photoelec- 
tric cross-sections given by Morrison & McCammon 
(1983). 



4.2. Energy Spectra Inside and Outside Dips 

We make use of joint fits of spectra from inside and 
outside a given dip to simultaneously constrain the 
parameters of the various spectral components. The 
total x 2 value obtained when the model is fit to multi- 
ple spectra is minimized by varying the model param- 
eters, but allowing only the absorption column den- 
sity of the bright component to vary among spectra. 
Since the disk blackbody and blackbody parameters 



typically evolve on timescales of hours (see Paper III ) 



it is important in this analysis to use a spectrum from 
immediately before (or after) the dip. 

The spectral fits were carried out on background- 
subtracted PCA pulse-height spectra. A 1% system- 
atic error estimate was added in quadrature to the 
estimated statistical error (la) for each channel of 
the spectra to account for calibration uncertainties. 
Fits were carried out only on data from PCUs 0, 1, 
and 4 and for energy channels corresponding to 2.5- 
25 keV, since the response model (version 2.2.1) for 
the other two PCA detectors and for outside that en- 
ergy range has been d etermined to be less accurate 
( Rcmillard et al. 1998| ). Fit parameters reported are 
the average values for PCUs 0, 1, and 4. Errors are 
conservatively estimated as the entire range encom- 
passed by the 90% confidence intervals from each of 
the detectors. 

Time segments immediately prior to dips 1 and 2 
( "A" ) and during the lowest portions of dips 1 and 2 
( "B" ) were selected for spectral analysis (dashed lines 
in Figs. H and ||). The two segments from Dip 1 were 
each 96 s in duration, and those from Dip 2 were 304 s 
in duration. The pulse-height spectra from these two 
time segments are shown in Figure ^ along with best- 
fit models. The parameters for the bright and faint 
component of each best- fit model are listed in the first 
two lines of Table [l]. The distance to Cir X-l has been 



estimated to be about 6-10 kpc (Stewart et al. 1991 
Goss fc Mebold 1977] ), so we chose a value of 8 kpc 



in converting blackbody and disk blackbody normal- 
izations to radii. The column density (Nh^) of the 
faint component is given in the table. The column 
density (Nh ) of the bright component differs for 
spectra A and B of each fit and is thus given sepa- 
rately in Table |[ 
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The disk blackbody and blackbody parameters for 
both dips are similar to those obtained for non-dip 
spectra in Paper III (temperatures of 1.5 keV and 
2.2 keV respectively). The unabsorbed flux of the 
faint component is ~10% that of the bright compo- 
nent. The column density of the faint component is 
very low, and consistent with zero. This low value 
is inconsistent with estimates of the interstellar col- 
umn density (Nh = 1.8-2.4xl0 22 cm -2 ) measured 
wit h ASCA and ROSAT ( [Brandt et al. 1996| ; |Prcdchl 
& Schmitt 1995 ). Apparently, the spectral shape of 



the faint component is slightly softer than that of the 
bright component, giving the appearance of reduced 
interstellar absorption. In fact, the ASCA data from 
a dip egress also showed that the low-state spectrum 
was softer than the high-state spectrum, and Brandt 
et al. (1996) hypothesized that this may be due to 
X-rays scattered into the line of sight by interstel- 
lar dust. This scattering by dust is distinct from the 
scattering by circumstellar material that might be re- 
sponsible for the low-state spectrum itself. The use of 
different spectral models for the bright and faint com- 
ponents or the addition of another continuum compo- 
nent would not be productive due to the large number 
of parameters involved, so we allow the column den- 
sity of the faint component to compensate for small 
differences in spectral shape. 

The column density of the bright component out- 
side dips (spectra 1A and 2A; see Table ||) is about 
3 x 10 22 cm -2 , just slightly above the interstellar value. 
The column density during the low segments of dips 
(spectra IB and 2B) is extremely high, i.e., 176 xlO 22 cm" 
for spectrum IB and 306x 10 22 cm -2 for spectrum 2B. 
This high column density value is required to produce 
the observed reduction in flux at 20 keV and is suffi- 
cient to render the absorbed component totally neg- 
ligible at lower energies. Thus, the flux in spectra IB 
and 2B below 5 keV may be attributed entirely to the 
faint component. 

4.3. Iron Emission Line 

Pulse-height spectra from the bottom of dips (spec- 
tra IB and 2B in Fig. Q) show a prominent peak near 
6.5 keV, suggesting iron Ka emission. The residuals 
for fitted spectra from both inside and outside the 
dips are similar, and in particular show peaks near 
6.5 keV with similar absolute flux levels. This sug- 
gests that the line feature is associated with the faint 
component. If the faint component is visible primar- 
ily because of scattering, iron fluorescence would be 



expected to occur in the scattering medium (see dis- 
cussion section below). 

The fits were repeated including a Gaussian emis- 
sion line component at 6.4-6.6 keV. Allowing the line 
flux to differ between fits of spectra A and B of each 
dip confirms that the line flux is the same to within 
a factor of two or less, despite an order of magnitude 
difference in continuum flux. Thus we jointly fit the 
spectra of each dip using the same line parameters for 
both spectra (Fig. ||). The addition of the line does 
reduce the magnitude of the residuals near 6.5 keV, 
but does not improve the fit below ~5.5 keV. The re- 
sults of the joint fits, including the Gaussian line, are 
listed in the second group in Tables | and |. They 
are generally consistent with the model parameters 
from the previous fits within the 90% confidence lim- 
its. Table [l] also shows that a significant improvement 
in the reduced \ 2 value was achieved in both cases. 
We will continue to include a Gaussian line in all sub- 
sequent fits and list its fitted parameters in Table [jj 
Although the best-fitting centroid energy of the line 
is close to 6.6 keV, the relatively coarse energy res- 
olution of the PCA (-1 keV FWHM at 6 keV) is 
insufficient to confirm whether the line is from neu- 
tral or partially ionized iron or whether it is actually 
composed of multiple unresolved narrow lines. 

4.4. Spectral Evolution During Dip Transi- 
tions 

Through its large collecting area, the PCA provides 
good-quality spectra of bright sources such as Cir X-l 
every 16 s. This enables us to study the detailed evo- 
lution of the spectrum during the transitions between 
high and low flux states outside and inside dips. We 
selected four to six 16-s segments from the ingress of 
dips 1, 2 and 3 (indicated as circles on Figures [|, [| 
and ||) and jointly fit the 4-6 spectra of each dip with 
the model developed above. The 16-s pulse-height 
spectra from each of these dips are shown in Figures 
[| [U^, and [n], together with the model curves and 
residuals. 

The low-energy intensity in each case initially de- 
creases but then reaches a fixed level while the inten- 
sity at higher energy continues to decrease. A "step" 
near 7.1 keV is evident in spectra (most prominently 
in spectra lc, 2c, 3d) for which significant, but not 
total, absorption of the bright component occurs at 
that energy. This step cannot be accounted for solely 
by the emission line but is naturally fit by iron K- 
edge absorption intrinsic to neutral or nearly-neutral 
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absorbing material with a column density (assuming 
solar abundances) of Nu > fewxlO 23 (intermediate 
and lower spectra, Table ^). As noted above, the 
intensity at the bottom of the dip is significantly re- 
duced relative to outside the dip even at 20 keV. 

The fit parameters for the bright and faint com- 
ponents in each joint fit of the 16-s spectra (la-d, 
2a-d, and 3a-f) are given in Table |. The best-fitting 
parameters for the spectra during dips 1 and 2 are 
very similar to the values obtained above for the fits 
(with iron line included) of the longer-exposure pairs 
of spectra, 1A/B and 2A/B respectively, from the 
same dips. The unabsorbed flux of the faint com- 
ponent relative to the unabsorbed flux of bright com- 
ponent is 10-12% in the fits of spectra la-d and 2a-d 
and somewhat lower, 5.4%, for spectra 3a-f. The col- 
umn density for the faint component is again consis- 
tent with zero in all cases. The Gaussian line param- 
eters for these fits are listed in Table [| and indicate a 
line centered at 6.4-6.6 keV. The best-fitting line flux 
differs by a factor of two among the three dips but 
is consistent, within the errors, with no change. For 
each set of spectra, the variable absorption column 
density of the bright component is shown in Table ||. 
In each case the column density increases from a rel- 
atively low value, N H = (4-9) xlO 22 , at the start of 
a dip to very high values, Nh > 10 24 , at the lowest 
part of the dips. 

These fits show that the model developed above 
is quite successful in reproducing the evolution of 
the spectrum through the dips. This justifies our 
assumption of common spectral values in the joint 
fits, where only the absorption column density on the 
bright component varies between all spectra of a given 
dip. 

5. Discussion 

During absorption dips in Cir X-l, about 10% of 
the total 2.5-25 keV flux remains unabsorbed. We 
conclude that the absorbing material does not com- 
pletely cover the X-ray source, and suggest that it 
is likely that the material totally covers our line of 
sight to a compact X-ray source during deep dips 
while a more extended source of X-rays remains fully 
or nearly fully visible. Similar behavior is observed 
during dips in other X-ray binaries. For example, 
an unabsorbed spectral component is detected dur- 
ing absorption dips in the binary pulsar Her X-l and 
is attributed to scattering in an extended corona (see 



Stclzcr ct al. 1999| and references therein). Likewise, 
curved CD/HID tracks similar to those we observe in 
Cir X-l were observed in the black- hole candidates 
GRO J1655-40 and 4U 1630-47 and were successfully 
modeled by absorption plus a faint unabsorbed com- 
ponent (Kuulkers et al. 1998; Tomsick et al. 1998| ). 



Our results are generally consistent with the ASCA 
results of Brandt et al. (1996, see above) for Cir X- 
1, but with the high-statistics RXTE data, we are 
able to demonstrate the detailed evolution of the con- 
tinuum spectrum through selected high-to-low flux 
transitions (dips). The ASCA spectra showed an 
iron K edge due to absorption with column densi- 



ties near 10 cm 



At 6.4 keV in the low-state 



ASCA spectrum, an iron line was seen with a flux of 
(8.9^5 6) xl0~ 3 photons cm -2 s -1 , and in the high- 
state spectrum, an upper limit of 4.5xl0~ 3 pho- 
tons cm~ 2 s _1 was derived for such a line. This in- 
dicates that the equivalent width of the line is much 
smaller in the high state, due to the brighter con- 
tinuum, but does not prove that the flux of the line 
was significantly smaller. Our analysis shows that a 
prominent iron emission line is present outside dips at 
a similar flux level, (12 — 25) x 10~ 3 photons cm -2 s _1 
(see Table ||), as inside dips. The equivalent width is 
therefore much smaller outside dips, as low as 44 eV, 
than during the deepest dips, where the equivalent 
width is typically 410-940 eV. Since scattering by 
an extended cloud of material is probably responsible 
for the faint component, it is likely that the iron Ken 
emission is also produced in the scattering medium 
(see below). 

If the scattering medium was a spherically symmet- 
ric cloud around the central X-ray source, scattering 
would not significantly alter the total flux of the sys- 
tem but would simply enlarge the size of the emission 
region for some fraction of the flux. We observe that 
~10% of the flux is in the diffuse component, requir- 
ing a column density of almost 2xl0 23 cm~ 2 to pro- 
vide sufficient Thomson-scattering optical depth. The 
iron line associated with the scattered component in- 
dicates that the scattering medium is not pure hydro- 
gen. Therefore, if the scattering material is not highly 
ionized, the direct component should show evidence 
for photoelectric absorption with a similar column 
density. The observed column density outside dips 
is about 2xl0 22 cm -2 — an order of magnitude less 
than required in this simple model. Thus we can con- 
clude that the scattered component is not produced in 
a spherical cloud of low-ionization state matter which 
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totally surrounds and covers the central source. 

The column density derived from the strength of 
the faint component can be used to estimate the iron 
Ka fluorescent flux in the scattering cloud. Based on 
the spectrum we measure above the iron K-shell ion- 
ization energy of 7.1 keV, a cloud of column density 
2xl0 23 cm~ 2 will photoelectrically absorb a flux of 
about 0.12 photons cm~ 2 s _1 . The Ka fluorescence 



yield for iron is about 0.3 (Bambynek et al. 1972), so 
we expect a flux in the emission line of 36xl0~ 3 pho- 
tons cm~ 2 s _1 . This is actually somewhat higher, by 
a factor of 1.5-3, than the line fluxes we have mea- 
sured (Table ||); however, this simple model illustrates 
that the iron Ka emission in Cir X-l is naturally ex- 
plained by the strong scattered component revealed 
during dips. 

We do not expect to find QPOs, all faster than 
1 Hz in Cir X-l, during portions of dips when the 
direct component is almost entirely absorbed, unless 
the scattering occurs in a region significantly smaller 
than a light-second (roughly the size of the accretion 
disk, 10 10 -10 n cm, discussed below). Power spec- 
tra during most of this 2-day observation show weak 
HBOs near 30 Hz, broad NBOs at ~4 Hz, or VLFN. 
No QPOs are detected during the lowest portions of 
absorption dips; however, the upper limits on QPOs 
in these low count rate circumstances do not allow us 
to exclude QPOs at a level similar to that outside the 
dips. 

The fact that absorption dips in Cir X-l mainly 
occur within a day of phase zero and are then fol- 
lowed by significant flaring suggests that the dips are 
associated with the mass transfer process. Thus, ab- 
sorption might be due to the mass transfer stream 
itself, a bulge or "splash point" near the edge of an 
accretion disk due to collision of an incoming stream 
with the disk, or a blob of outer disk material pro- 
truding vertically into the line of sight. Although ab- 
sorption dips in Cir X-l show a wide range of profiles 
and durations, for a brief deep dip we can attempt to 
esti mate the size and density of an absorber near the 



lib 



out er edge of an accretion disk, 



We first consider the simple case of blobs in Ke- 
plerian motion near the outer edge of a disk. Other 
than the 16.55 d orbital period of Cir X-l, the bi- 
nary parameters, including a probable high eccen- 
tricity, are unknown. However, the Keplerian ve- 
locity in the outer edge of a disk does not depend 
strongly on the binary parameters, so we adopt val- 
ues of Mns ~ Mdonor ~ 1.4 M© and eccentricity ~ 



0.8. An outer disk radius of R ~ 10 10 -10 11 cm, giv- 
ing a Keplerian velocity of order 1000 km s _1 , would 
then be about 5-50% of the radius of the instanta- 
neous critical potential lobe of the neutron star near 
periastron. A typical brief deep dip such as Dip 1 
lasts of order 100 s, implying an azimuthal dimension 
of ~T0 10 cm, or a significant fraction of the accre- 
tion disk size. If the radial and azimuthal dimen- 
sions of the blob are similar, a column density of 
10 24 atoms cm -2 requires a number density in the 
blob of n ~T0 14 atoms cm~ 3 . For typical luminosities 
of L x ^3xl0 38 erg s _1 cm~ 2 (determined from the 
unabsorbed flux, assuming a distance to Cir X-l of 
8 kpc), the ionization parameter £ = L x /nR 2 for a 
blob of this density would be ~10 3 -10 4 , which would 
in turn imply that the blob of gas wou ld be nearly 
fully ionized ( Kallman fc McCray 198*2] ). However, 
such a highly ionized absorber is inconsistent with 
the strong Fe K edge near 7 keV observed during dips 
in Cir X-l. The large size of the blobs in this picture 
relative to the disk size might indicate that the blobs 
are actually azimuthally extended or that dips, par- 
ticularly the longer dips, are due to multiple blobs. 
The Keplerian period of blobs in this picture is of or- 
der 10 3 s; however, no recurrence of dips on this time 
scale is apparent in the light curves, suggesting that 
such blobs may be short-lived. 

Alternately, we consider absorbers that are co- 
rotating with the binary orbit, e.g, a stationary "splash 
point" on the outer edge of the disk. Using the same 
simple assumptions about the binary parameters and 
disk size as above, we find an azimuthal velocity of 
order 10 km s _1 for such an absorber. This slower 
velocity would require a much smaller absorber, with 
dimensions ^10 8 cm, compared to a blob in a Keple- 
rian orbit. The corresponding number density in the 
absorber would therefore be ~10 16 atoms cm -3 . The 
ionization parameter would be lower for these denser 
blobs, £ ~ 10-100, and the absorbing gas would not 
be fully ionized, consistent with the photoelectric ab- 



sorpt ion observed during dips ( Kallman fc McCray 



1982J). In this picture, the fact that Cir X-l actu- 



ally exhibits many dips during a couple days of each 
binary cycle suggests a stationary (in the rotating 
frame) structure in the accretion disk (i.e., a station- 
ary wave) or else complex structure (e.g., shocks or 
turbulence) in the accretion stream. 

The chaotic nature of the PCA light curves during 
dips (Figs. 0-|J) indicate that we are observing fine 
structure (temporal or geometrical) of the obscuring 
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material. The true nature of the absorbing material 
is likely to be as complex as the light curves them- 
selves, reflecting the turbulent conditions associated 
with mass transfer in an eccentric binary. In Cir X-l 
the accretion disk probably experiences a high de- 
gree of disturbance due to a high eccentricity with 
enhanced mass transfer occurring near periastron. 

We acknowledge the RXTE teams at MIT & GSFC 
for their support. In particular we thank E. Morgan, 
R. Remillard, W. Cui, and D. Chakrabarty for use- 
ful discussions. Support for this work was provided 
through NASA Contract NAS5-30612. 
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Fig. 1. — Light curves in four energy channels and two hardness ratios for PCA observations of Cir X-l from 
1996 September 20-22, covering a two-day period around phase zero (0 = 0). Three dips have been identified for 
further study. Each point represents 16 s of background-subtracted data from all five PCA detectors. Ratios of 
the intensities produce soft (I [4.8-6.3 keV] / I [2.0-4.8 keV]) and hard (I [13-18 kcV] / I [6.3-13 keV]) hardness 
ratios. The intensity levels of the segment at day ~347.2 (after Dip 1) are close to the level in each band during 
the quiescent phases of the orbit. 10 
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Fig. 2. — Expanded 2400-s segment (1-s bins) of the light curves and hardness ratios in Fig. [I] including Dip 1. 
Time zero corresponds to day 347.043 in Fig. [l| Energy spectra were extracted from two 96-s time segments (A 
and B), indicated by dotted vertical lines, and four 16-s segments (a-d) indicated by circles. 
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Fig. 3. — Expanded 2400-s segment (1-s bins) of the light curves and hardness ratios in Fig. [I] including Dip 2. 
Time zero corresponds to day 347.272 in Fig. |l|. Energy spectra were extracted from two 304-s time segments (A 
and B) and four 16-s segments (a-d). 
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Fig. 4. — Expanded 2400-s segment (1-s bins) of the light curves and hardness ratios in Fig. [I] including Dip 3. 
Time zero corresponds to day 347.344 in Fig. |l|. Energy spectra were extracted from six 16-s time segments (a-f). 
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intensity gradually transitioned from the non-dip baseline to the bottom of a dip. Circled points labeled 'a-f 
correspond to the six 16-s time segments identified in Fig. || and used to construct energy spectra. Intensity 
is I [2.0-18 keV], and the hardness ratios are defined as soft color: I [4.8-6.3 keV] / I [2.0-4.8 keV], hard color: 
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background-subtracted data from all five PCA detectors. 
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. 9. — Spectral fits for four 16-s segments (a-d) during the decline into Dip 1. Only the column density of the 
jht component varies between the four jointly-fit curves. The data shown are from PCU only. 
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Spectral fits for four 16-s segments (a-d) during the decline into Dip 2. See the caption of Fig. ^. 
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Table 1 

Joint fit parameters of the bright and faint spectral components during dips. 



Disk Blackbody Blackbody Faint Component 

Joint Iron kT in R m cos 1 / 2 B c kT R d /10 22 f 

fit a Line b (keV) (km) (kcV) (km) / c (cm" 2 ) xl 



1A-B 
2A-B 



no 



1.51 
1.51 



+0.10 
-0.17 
+0.04 
-0.06 



16.91 



+3.96 
2.20 



2.20 
2.32 



-0.62 
-0.42 

-0.16 
-0.15 



3.77 
3.47 



+3.33 
2.55 
+0.86 
0.76 



Q 107 +0.005 
u - lu '-0.006 



0.103 



+0.003 
0.003 



,+0.07 
J -0.00 

0.03+nn? 7.58-7.84 



1A-B 
2A-B 



yes 
yes 



1.23 
1.35 



+0.22 
-0.15 
-0.06 
-0.05 



25.32 



-8.56 
-7.03 



22 4^+ 1 ' 83 



1.79 
2.10 



-0.33 
-0.16 
+0.12 
0.09 



7.87 
5.23 



+3.34 
4.14 
+0.92 
0.93 



096+ ' 007 
u.u»o_ 006 

092+ 002 



-,+0.10 



0.00 



+0.10 
0.00 



1.87-2.15 



la- 


-d 


yes 


1 - iS -0.16 


27 04+ 9 39 
Z(.u*_ 6 33 


i 79 +0.22 
1 -' z -0.12 


o 99 +2.69 
°' zz -3.03 


Q 10 o+0.009 


00+ ' 44 
u - uu -o.oo 


0.99- 


-1.20 


2a 


d 


yes 


1 26+ 014 
1 - ZD -0.14 


2 c i 2 +7.66 


2.08±°;f 5 


5 46+ 1 0S 
°- TO -1.08 


n7 +0.013 
u - il( -0.014 


15+ - 70 
U - i0 -0.15 


0.80- 


-1.12 


3a 


-f 


yes 


1 07 +o.08 

1 - u '-0.09 


42 66+ 8 ' 34 


1 74+0.09 

1 -' y: -0.09 


9 23+ 1 ' 41 
a,zo -1.32 


054+ 005 


00+ 14 
u - uu -o.oo 


1.27- 


-1.34 



a Joint fits of spectra from outside (A) and inside (B) dips 1 and 2 and for the four to six 16-s spectra (a-f) 
during dips 1, 2, and 3 (using the model given in equation |lj). |-Pnly the column density (N^) of the bright 
component is allowed to vary between spectra (see Tables U and W). Errors quoted are 90% confidence limits for 
a single parameter (Ax 2 = 2.7). 

b Gaussian emission line at 6.4-6.6 keV was included in the final five fits (see Table ^|). 

c Inner radius of the accretion disk (times cos 1 / 2 8, where 6 is the angle between the normal to the disk and the 
line of sight) for a distance of 8 kpc. 

d Blackbody radius for a distance of 8 kpc. 

e Ratio of the unabsorbed flux of the faint component to the unabsorbed flux of the bright component. 

f Absorption column density of the faint component. See Table ^ for the variable absorption of the bright 
component. 

g x 2 /^°/i where dof = the number of spectral bins (52-54 per spectrum) minus the number of fit parameters. 
The range of values represents the fits for PCUs 0, 1, & 4. 



21 



Table 2 

Absorption column density of the bright component outside (A) and inside (B) Dips 1 and 2. 



A^/10 22 



Spectrum (cm 2 ) Spectrum (cm 2 ) 





No Emission Line 




1A 


2.83±g;f| 2A 


o cn+0.29 
°- ou -0.34 


IB 


176^13 2B 


306t 2 g 




With Emission Line a 




1A 


3.81±g:S 2A 


4 oq+0.26 


IB 


184+^ 2B 


283±^ 



a A Gaussian emission line at 6.4-6.6 keV was in- 
cluded, with the same parameters outside and inside 
the dips (see Table ||) . 



Table 3 

Gaussian emission line parameters for the joint fits of dip spectra. 



Joint 


E 




Line flux 


Fit 


(keV) 


(keV) 


(photons cm' 2 s _1 ) 


1A-B 
2A-B 


e^ 2 
6.58ir 6 7 


0.181°;?! 
0.29l°;i° 


o.oisi™ ,! 

O.Oieiffi 


la-d 
2a-d 
3a-f 


6.45+8:2 

6.46t°; 23 7 


u - 1,3 -0.13 

u - 4 ' -0.32 

QQ+0-32 
U.O3_ 3g 


U.U10_Q QQ2 

u.uiz_ 003 



Gaussian width. 
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Table 4 

Effective hydrogen column density responsible for the variable absorption of the bright 
component during 16-s spectra from dips 1, 2, and 3. 



JV^/IO 22 AT«/10 22 AT^/10 22 

Spectrum (cm -2 ) Spectrum (cm -2 ) Spectrum (cm -2 ) 



la 


4 1fi+ 110 


2a 


8.84tl;jg 


3a 


o O1+109 
°-' 31 -0.96 


lb 


25.8±1; 6 3 


2b 


25.8t\i 


3b 


15.7±i;? 


lc 


43.6^;I 


2c 


64.0tH 


3c 


33.2±i| 


Id 


129±| 


2d 


242±f 2 


3d 


66.7±1;| 



3e 1271? 
3f >2xl0 4 
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